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ABSTRACT 


A mathematical model of the annoyance created at an airport by air- 
craft operations is developed, The model Incorporates population distri- 
bution considerations around an airport and the annoyance caused by air- 
craft noise. The objective function of this model corresponds to seeking 
to minimize total population annoyance created by all aircraft operations 
in a 24-hour period, Several factors are included in this model as con- 
straint relationships. Aircraft operations by type and time period are 
upper bounded. Demand for flight services is incorporated by including 
lowei" bounds on the number of operations by type of aircraft, runway used 
and time period. Also upper bounds on the number of operations for each 
runway are included. The mathematical model as formulated is recognized 
as corresponding to a nonlinear integer mathematical programming problem. 
The solution technique selected makes use of a successive linear 
approximation optimization algorithm. An especially attractive feature 
of this solution algorithm is that it is capable of obtaining solutions 
to large problems. For example, it would be feasible to attempt the 
solution of problems involving several thousand variables and 500 plus 
linear constraints. This suggested solution algorithm was Implemented on 
a computer and computational results obtained for example problems. 


I , INTRODUCTION 


With the advent in the late 1950 's of the jet~engined commerlcal air- 
craft, airports began facing an increasing noise pollution problem. As 
aircraft si.ie and range increased there was a corresponding increase in 
the number of commercial routes over which these carriers operated. More 
routes dictated more aircraft, so existing airports were required to han^^le 
an increased volume of air traffic each year. New airports were constructed 
to meet the ever increasing demand for air transport, and the older, less 
scientifically designed airports continued to handle higher and higher 
volumes of traffic. 

Airports were designed to handle not only Increased volumes of air 
traffic but also larger aircraft. What the airports, up to this point, 
have not been able to satisfactorily handle is the associated noise pollu- 
tion. As the noise pollution has increased, so has the opposition to the 
noise. This opposition has become well organized and has addressed not 
only noise levels within airport boundaries, but noise pollution in the 
acreage surrounding airports as well. Noise, or more specifically the rt 
ductlon of noise around airports has become a critical economic as well as 
emotional issue. 

As with most of today's complex issues there is not an easy solution 
to the problem of aircraft noise. In fact, there is not even a consensus 
of opinion as to what constitutes the best solution. Individuals who are 
annoyed by the noise argue that the noise is Interferring with their lives 
and affecting their property values; they implore the airports and airlines 
that fly into them to reduce the noise to an acceptable level immediately. 
The airlines argue that any short-term solution to the problem would be a 



disaster to them economically and urge understanding until they can econom- 
ically phase out their older, noisier jets and replace them with newer, 
quieter models, Airport managers say they are caught in the middle, sym- 
pathizing with the individuals who are annoyed by the noise, but powerless 
to individually effect a change on the airlines without strained relations 
and hard feelings. 

The goal of this research is to determine the minimum noise impact 
assignment - scheduling of aircraft to the existing arrival and departure 
trajectories for any airport of concern. To support the modelling effort 
for the minimum noise formulation a considerable amount of background 
material is required. The next two sections contain this related material. 
The problem of interest is formulated via a mathematical model in section IV. 
The objective function of this model corresponds to seeking to minimize the 
total population annoyance metric as a function of aircraft operations in 
a 24-hour period. Factors such as aircraft operations, demand and time 
restrictions are Included as constraint relationships. Section V presents 
a solution technique that makes use of gradient Information and successive 
linear programming to obtain approximate noise minimal solutions. An 
attractive feature of this solution algorithm is that one can efficiently 
use it to produce solutions and sensitivity analyses for very large pro- 
blems. Computational results are then presented for an application air- 
port. These results are especially encouraging when compared with current 
operating scenarios. That is, they indicate a substantial reduction in 
airport community annoyance may be achieved by merely changing the assign- 
ments of servicing aircraft to the established arrival and departure 
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II, BACKGROUND INFOW-IATION ON NOISE 


II. 1 Sound Measurement: 

The problem associated with the measuring of sound around airports is 
one of developing a rating scheme that will quantify the sound as acceptable 
or unacceptable to the human ear. The basic instrument used to measure 
sound is the sound level meter. If this sound level meter is used to mea- 
sure sound regardless of frequency, called the Overall Sound Pressure Level, 
a measure called the C-Weighted Sound Level is derived, This meter permits 
one to obtain a reading that will be satisfactory in determining the sound 
acceptability and in rank-ordering various noises, Although the C-Weighted 
Sound Level is a measure of sound it does not correspond very well to the 
way individuals judge noises [12], 

A measure of sound that more accurately corresponds to the way people 
judge noises is ‘■he A-Weighted Sound Level (L^) , or A-Level, This measure 
is expressed in decibels dB(A) and is the weighted sum of all of the compo- 
nents of the noise. This measure has been found to correlate very well to 
an individual’s subjective judgment of the annoyance of many types of noise. 
Table 1 gives representative values of A-Weighted Sound [12J . 

11,2 Human Reaction to Noise 

Much more difficult than measuring the sound Is quantifying the human 
reaction to a measured noise level. Noise can have one ot three effects on 
people: (1) Subjective (annoyance, nuisance, dissatisfaction, disturbance, 

etc.); (2) Behavioral (interference with an on-going activity); and (3) Physio- 
logical [12] . Table 2 depicts the noise levels that will most likely inter- 
fere with specific activities. 

Even being able to measure noise satisfactorily and then knowing what 


TABLE 1 


Decibel Levels 


Measure dB(A) 

Perceived Loudness 

0-10 

Tlareshold of hearing 

10-40 

Barely audible 

40-60 

Quiet 

60-80 

Moderately loud 

80-110 

Very loud 

110-130 

Uncomfortably loud 


Source ; U.S. Department of Housing and Urban 
Development. "Noise Assessment 
Guidelines Technical Background." 
Report //TE/NA 172, Washington, D.C. 


TABLE 2 

Noise Interference 


Activity 

Average, outdoor 
noise levels 

Interference with Radio & 
T.V. listening with the 
windows open 

45 to 50 dB(A) 

Interference with Radio & 
T.V. listening with the 
windows closed 

55 to 60 dB(A) 

Sleep interference 

40 to 45 dB(A) 

Acceptable living 

Less than 

environment 

80 to 85 dB(A) 


Source : U.S. Department of Housing and Urban 

Development. "Noise Assessment Guidelines 
Technical Background." Report //TE/NA 172. 
Washington, D.C. 
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levels of that noise will Interfere with everyday functions is not totally 
sufficient. Studies have shown that approximately lOZ of the population 
are so sensitive to noise that any noise not of their own making is objec- 
tionable to them. Conversely, it has been estimated that 25% of the popula- 
tion is virtually Insensitive to noise and will not complain even in very 
severe noise environments [12] . 

Individual reaction to noise will depend on various items including 
previous noise exposures, psychological attitudes, socio-economic status 
and the nature of the activity that is intruded upon by the nolsev hue 
to the disparity in these factors, any derived scale will be subject to 
inaccuracies and question. It is more important to look on a rating scale 
(or metric) as an attempt to describe, in physical terms, the nature and 
magnitude of the total noise as it affects groups of individuals. The 
scale is helpful in determining the level of noise stimulus that cannot be 
exceeded without rendering the environment unacceptable for living for most 
people [12],, Any scale will attempt to account for the context in which 
the noise stimulus is experienced and can well Introduce various adjustments 
for characteristics of the noise and the situation on which the noise 
intruded, 

II. 3 Noise Measurement and Assessment 

As has been discussed, a mere fluctuation in the needle of a Sound 
Level Meter cannot indicate the level of annoyance that particular sound 
will create in an individual. To quantify annoyance, noise exposure must 
be introduced. Noise exposure is defined as "the whole time-varying pat- 
tern of the sound level rather than some single level, such as the average 
value" [12] . Results of experiments show that a steady noise is more 
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acceptable to people than a noise of the same average level that fluctuates 
erratically, In most cases, the greater the fluctuation, the greater the 
annoyance. Several single-event noise measurement scales will be introduced 
and their relationship to the measurement and quantification of multiple- 
event scales, that can be used to assess the noise environment around an 
airport, will be specified. 

11. 4 Perceived Noise Level (PNL) 

The PNL may be defined as "a quantity that is calculated from measured 
noise levels and adjusted by weighting more heavily those frequencies that 
are more annoying to the listener" [7], 

PNL « 1.02 L^ + 11.5 (1) 

Where L^ represents the noise level in dB(A). The PNL is measured in PNdB. 
This measure is based on the Individual judgments of "equal annoyance for 
bands of sound one-third octave wide during an aircraft flyover" [11] , 

11. 5 Effective Perceived Noise Level (EPNL) 

This scale, expressed in EPNdB is defined as a unit of perceived noise 
that "takes into account the actual sound energy received by a listener, 
the ears' response to that sound evergy, the added annoyance of any pure 
tones or 'screeches’ in the noise, and the duration of the noise" [13],. 

The EPNL has been adopted by the Federal Aviation Administration (FAA) as 
a measure of the noise emission level of individual jet aircraft. 

11. 6 Noise Exposure Level 

Another single event noise metric is the Noise Exposure Level (NEL) . 

The NEL may be defined as the summation of the time varying sound level 



over the time span the sound level is within 30 dU of the maximum. In 
equation form 

t2 (L (t)/10) 

NEL ° 10 log Q / 10 ^ dt (2) 

where tj^ and correspond to sound levels of 30 dB less than the maximum 
sound, corresponds to the A-level sound at time t. 

11,7 Noise Exposure Forecasts 

A method of measuring noise from multiple events around airports is 
called the Noise Exposure Forecast (NEF) . The NEF may he defined as "the 
computed summation over a 24-hour period based on Effective Perceived Noise 
Level, the number of daytime (7 a.m. to 10 p.m.) 'ircraft noise events, 
and the number of nighttime (10 p.m. to 7 a.m.) aircraft noise events" [6]. 
This forecast "provides a measure of the total aircraft-generated noise 
received at locations near an airport during a typical 24-hour period" [13]. 
The NEF value may be represented mathematically as: 


NEF, . 


= (EPNL) , + 10 log (N + 17 ) ~ 88 

^ ij ”ij 


(3) 


where the terms are defined as: 

EPNL = effective perceived noise level for the particular 
aircraft at the given point on the ground, in 
decibels (A) 

1 = aircraft class 

j = aircraft flight path 

N^ = number of daytime events of the particular aircraft 

ij 


= number of nighttime events of the particular aircraft 


*• 8 •“ 

Once these ^.ndividual NEF values are computed, the total NEF value 
of a ground location is computed by the equation: 


NEF « 10 log Z S antilog (NEF, J/10 (4) 

i j 

Table 3 represents typical NEF values and the level of public complaint 
that such values will evoke. To scale the NEF values, one may think of the 
reduction of one NEF unit as the equivalent of the reduction of approximately 
2% of the number of people who will be highly annoyed by the noise [13] . 

Once the NEF values of locations surrounding airports are computed, 
locations of equal NEF may be joined together to give NEF isopleths. It 
should be noted that the Federal Aviation Administration (FAA) has accepted 
it for land use planning around commercial jet transport airports [12], 

Along with the FAA, the U,S. Department of Housing and Development (HUD) 
has published the guidelines in Table 4 as a site screening device for 
residential housing. 


11.8 Day-Night Level 

Day-night level (L^^) was developed as a single-number measure of 
community noise exposure. It is defined as the average A-Welghted noise 
level integrated over a 24-bour period. Appropriate we: ghtlngs are applied 
for the noise levels occurring in the daytime and nighttime periods [3]. 

It is stated in [15] that is "the primary measure for describing noise 
in an environmental impact statement". For discrete sampling of A-Weighted 
sound level for a 24-hour time period, may be formulated mathematically 


as. 


n 


Ldn = 10 log [• 


Z w, antilog (L. ,/10) 

i— 1 ^ AjX 


n 


] 


(5) 


m \ 


^ xM 

TABLE 3 

NEP Interpretation 


NEF Value 


Interpretation 


Less than NEF 30 


NEF 30 to NEF AO 


More than NEF 40 


No complaint expected, noise 
may interfere with community 
activities. 

Individual may complain, 
group action possible. 

Repeated vigorous complaints 
expected, group action 
probable. 


Source ; U.S. Department of Transportation. "Aviation 
Noise Abatement Policy," Washington, D.C., 
November 18, 1976. 


TABLE 4 

HUD NEF Interpretation 


NEF Level 


Residential Site Category 


Less than NEF 30 

Normally acceptable for any type 



of construction. 

Between 

NEF 30 and 

Normally unacceptable for single- 

NEF 40 


unit residential construction; 
acceptable for multi-unit 
construction with soundproofing. 

Greater 

than NEF 40 

Unacceptable for practically all 
types of residential construction 
area restricted to agricultural, 
outdoor recreational or industrial 
uses . 


Source ; Cawthorn, Jimmy M. and Brown, Christine G. "Effect of 
Advanced Aircraft Noise Reduction Technology on the 
1980 Projected Noise Environment Around Patrick Henry 
Airport." NASA Technical Memorandum. Langley Research 
Center. Hampton, Virginia, 1974. 
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where 

« time of day weighting for sample i 
La i “ A-level for sample i 

n ° number of samples of La in a 24-hour period. 


II. 9 Noise Exposure 

A set of computer programs has been generated for the U.S, Department 
of Transportation for the computation of noise exposure values due to air- 
craft operations around airports. The collection of programs is called the 
Integrated Noise Model (INM) , The INM system is available from the Federal 
Aviation Administration (FAA) , Washington, D. C., at a nominal charge. The 
INM will compute noise exposure values for the following noise metrics; 
Noise Exposure Forecast (NEF) ; Day-Night Sound Level (L^^^) ; Community Noise 

Equivalent; Level (CNEL) ; Equivalent Sound Level (L ) ; Aircraft Sound 

eq 

Description System (ASDS); and an experimental metric called Dally Overall 
Sound Exposure (DOSE). The use and description of INM is presented in 
references [1] and [4]. The calculations of NEF and L^^ noise exposures 
can be represented by the following generalized equation; 

n (EL /lO) 

NE = 10 log { Z (aD, + bE + cN.) 10 ^ } - A (6) 

i=l ^ ^ ^ 


where, 

represents noise exposure (either NEF or L, ) 

dn 

a, b, c day, evening and night weighting factors (for 
NEF and L^^, c “ 16.7 and c = 10 respectively 
while a = b = 1 for both metrics) 

D, E, N actual number of day, evening and night 
(respectively) flight operations 
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EL single event exposure level (i,e,, EPNL for 

NEF, or NEL for L^^^) 

A 88.0 for NEF or A 9. 4 for 

Equation (6) is a concise representation of the noise exposure calculation 
process. This calculation, briefly, consists of determining the relative 
position of an aircraft to a point of interest and the physical components 
computed (i.e., thrust end velocity). Then, the single event level is 
found and the weighted number of identical operotions of the type being 
considered are factored into the computations. Finally, a cumulative sum 
of noise exposure comprising the contributions from all distinctly different 
kinds of operations from all the aircraft flights on all the ground tracks 
is computed yielding the final total exposure, 

One 3.ast metric will be presented, This metric not only uses a mul- 
fjple noise metric but simultaneously weights the impact with popu- 

lation figures. This metric is then used in formulating an objective in 
the mathematical representation of the airport noise problem (section IV). 

11,10 Noise Impact Index 

The noise impact index (Nil) , may be used for comparing the relative 
impact of one noise environment with that of another. "It is defined as 
the sound level weighted population divided by the total population under 
consideration" [15]. The formula for this index is 


^Total (7) 

where 

LWP =* sound level weighted population 

Pm “ total population under consideration. 

Total 






The sound level populaclou represents the integrated effect of given noise 
environments on a particular population, LWP is represented mathematically 


LWP = / i?a<,„)-wa^„)<KL,„) 


wnere 


P(L^^) " population distribution function 

W(L^n^ ” day’-nlght average sound level weighting function 

d(L^n) “ differential change in day-night average sound level, 


An effective approximation for the calculation of the Noise Number Index is 


Nil « 


k 


where 


“ population in area k 


W((Lj^)j^) “ day-night average sound level weighting function 


^^dn^k 


day-night level for area k. 


An example of W(L ) would be the sound level weighting function for overall 
dn 

impact analysis described in’ [15] . The analytic expression for this function 


W(^dn> = 


, 0.103L 

[3.364 X lO”*"] [10 

0 . 03L , ' 7 0.08L, 

[0.2] [10 ^"^] + [1.43 X 10 ][10 ‘^^] 
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III. POTENTIAL STRATEGIES FOR REDUCING AIRPORT NOISE POLLUTION 

Just as with the building of airports and the establishing of federal 
airways, considerable thought needs to be given to the long terra effects of 
noise on populations living in proximity to major airports. The effect of 
aircraft noise pollution is not confined just to residents in the immediate 
area of the airport, as many complaints about noise come from people living 
at some distance from the airport. These people usually reside near either 
the approach paths into an airport or the take-off paths out of an airport. 
Proposed strategies for reducing airport associated noise usually consist 
of methods that airport operators may use to decrease the size of the NEF 40 
and NEF 30 contours surrounding their airports. Recently, there has been a 
substantial amount of work being done in the area of engine modification as 
a way of reducing aircraft noise. 

Flight procedures tend to be effective in reducing noise pollution 
but are controversial due to the safety aspects. Changes in flight proce- 
dures can affect either take-off or landing with the take-off involving 
primarily jet noise and the approach involving primarily machinery noise [7], 
"In spite of the fact that much higher engine thrust is required for take- 
offs than landings, landing noise is frequently more annoying to the ear 
because of dominant fan noise" [11] . Landing approaches also tend to be 
less steep than take-offs, so a greacer land area is exposed to this low 
latitude noise for a longer period of time [5]. 

Presently most aircraft use a 3° approach angle (one segment), or 
glide slope as it is called, which results in the aircraft being in its 
final landing configuration (flaps down, landing gear down) , stabilized in 
speed and power at a height of not less than 1000 feet above the ground. 
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A two segment approach calls for a 6° glide slope from 3500 feet to 1000 
feet. At 1000 feet a transition is made from a 6" glide slope to a 3“ glide 
slope, V7ith the transition complete at 700 feet above the ground [8] . The 
tv’O segment approach would Involve additional avionics in the aircraft and 
additional navigational aids on the ground, but would not require modifica- 
tion of any aircraft. Studies have stated that for the noisier aircraft in 
the current fleets, the two segment approach will reduce the NEF 30 area by 
26%, but for the newer quieter aircraft the two segment approach would have 
little effect on noise reduction [2]. The airlines and Air Traffic Control 
(ATG) have objected to the two segment approach as a safety hazard and have 
Instead devised a low drag/low power approach. This solution still uses a 
3° glide slope with interception at 3000 feet as before, but with an inter- 
mediate flap setting instead of full flaps. Landing gear is delayed until 
700 feet [8]. 

The most accepted method of noise abatement on take-off is to take-off 
under full power and climb at the steepest possible angle to gain height 
before flying over populated ares':;. Once the height is achieved to conform 
to safety procedures the engines may be throttled back during a flyover of 
populated areas. If this idea of a steep climbout is coupled with holding 
a constant speed instead of acceleration during climbout, a reduction of 
as much as .6 dB(A) may be realized [10]. Another method that may be used 
in some areas during climbout is the execution of turns in a direction aways 
from populated areas. This also assists in decreasing noise exposure and 
annoyance. One technique that may be used in good weather, and is presently 
used at Washington’s National Airport, is the concentration of take-offs to 
strictly defined corridors. These corridors correspond to the Potomac River. 
In this method fewer people are inconvenienced, but the leyel of annoyance 
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within these corridors is increased [5], If such corridors can be estab- 
lished over unpopulated, or favorably zoned land, then they can have a 
significant effect by decreasing the annoyance due to airplane noise. 

Several parameters that affect the airport noise environment directly 

are, 

(1) runway usage and trajectory selection for arriving and 
departing aircraft, 

(2) total number of daily operations, 

(3) operations by aircraft type and time period limitations. 

There is little that can be done about runway orientation once an air- 
port is operational. Even if an airport is in the design stages, safety 
considerations will dictate that runway orientation be in consonance with 
the prevailing winds in the area, Prevailing winds will also affect the 
utilization of runways, dictating on many occasions which runways may be 
used. When wind is not a factor, however, the option of runway and trajec- 
tory assignment may be exercised, 

Another option that may be exercised in attempting to reduce noise 
around airports is to limit the number of daily flights of certain types of 
aircraft. This is a logical step since the NEF contours that quantify annoy- 
ance are calculated based on the fleet mix as well as the number of flights 
into and out of an airport within a 24-hour period. 

Airport authorities cannot dictate the types of aircraft that utilize 
their airports, but they are aware that various types of aircraft have dif- 
ferent effects on the noise pollution around their airports. Initially, 
without limiting the tiumber of noisy aircraft, the managers may limit the 

time of day that a noisy aircraft may land. If the NEF or L, equations are 

dn 

examined, one will find that nighttime aircraft noise events are weighted 
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much greater than daytime aircraft noise events. If the number of noisy 
aircraft that are allowed to land at night are reduced, then the overall 
noise impact may be reduced, 

One form of time period limitations was introduced by National Airport 
in Washington, D. C. as a part of its 1978 Draft Environmental Impact State- 
ment [14]. Between the hours of 7 a.m. and 9 p.m., a maximum of 40 schedules 
per hour (take-off and landing) were proposed. From 9 p.m. until 9:30 p.m., 
20 schedules were proposed. After 9:30 p„m., the number and type of air- 
craft that may depart or land is severely limited. Merely examining the NEP 
or equations will Indicate that this action may reduce the noise exposure 
as measured by the response indices. 

It is these last three parameters with which this research is primarily 
concerned. Our proposed objective of minimizing annoyance from aircraft 
noise will be sought by selecting the optltpal trajectory selection and the 
optimal assignment of aircraft, by time period, servicing the airport of 
interest. This optimization will be conducted subject to constraints on 
demand for service, allowable number of operations, aircraft availability 
and geographical area allowable noise levels. The mathematical formulation 
of the problem of interest follows. 

IV. AIRPORT NOISE MINIMIZATION MODEL 

Assumptions that are made in developing the model for any airport of 
interest are: 

a. The model is developed primarily for application to commercial 
airline traffic. 

b. Approaching aircraft follow one of a group of fixed inbound 
trajectories. Likewise departing aircraft follow one of a 
group of fixed outbound trajectories. 
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The proposed model will be formed to Include the following control 
parameters ; 

a. The number of aircraft of each type assigned to each specific 
trajectory in any designated time period for take’-off, 

b. The number of aircraft of each type assigned to each specific 
trajectory in any designated time period for landing, 

The following symbols will be used in the formulation of the mathe-^ 
matical model: 

VARIABLES 

X 

Ikj/l number of departures of type i aircraft with stage 

length k utilizing trajectory j during period i 
Y 

IjJl = number of arrivals of type 1 aircraft utilizing 
trajectory j during period H 

CONSTANTS 

A “ area designation 

i = type of aircraft designation 

j = ground track designation 

k = stage length designation 

H != time period designation 

P = total affected population 

R = runway number 

NA “ number of areas 

ND = number of runways for departure aircraft 

NI = number of types of aircraft 

NJ’ - number of different ground tracks 

NK = number of stage lengths 


NL 


number of time periods 


NR 

NV 


N 


N. 


•A 


R 


j 


"a 

NX 


U 



tsdp 


’=> number of runways 

» number of runways for arriving aircraft 
“ limitation on number of operations that use ground 
track d 

« limitation on number of operations in time period p 
« limitation on number of operations for runway r 
>= limitation on number of operations with stage length s 
= limitation on number of airport operations per day 
» population in area A 

= set of ground tracks associated with runway R 
= community response index critical value for area A 
= number of “ype i aircraft available for take-offs 
during time period i 

- number of type i aircraft available for landing 
during time period i 

- number of type t aircraft with stage length s with 
ground track d' that are required during time period p 


NOISE PARAMETERS 


EPNL 


EPNL 


ijA 


EPNL 


IkJA 


L 


A 


effective perceived noise level 

effective perceived noise level for arriving aircraft 
of type i corresponding to ground track j experienced 
in area A 

effective perceived noise level for departing aircraft 
of type i with stage length k corresponding to ground 
track j experienced in area A 
A-Weighted sound pressure level 


° day-night level 
“ equivalent sound level 
NEF “ noise exposure forecast 

RI^ o coiranunity response index for area A 

W(RIy) » weighting factor as a function of community response 

index for area A 

IV. 1 Ob;jective Function 

The objective function will be developed using population Information 
and weighting factors discussed in Section II. For each of the grid areas 
shown in Figure 1, there is an associated population that is assumed to 
be evenly distributed throughout the grid. If the populations of all the 
grid areas are summed, then the resulting P will be the entire affected 
population around an airport. The fraction of the population affected in 
any grid area A is: 

P^/P (11) 

As is discussed in Section II. 9, the selected community response index 

for grid A, » will be utilized in conjunction with an associated weight- 

A 

ing factor for objective function formulation. The selected objective is to 
minimize the sound level weighted population divided by the total population 
under consideration. Mathematically, this is represented as: 

NA 

minimize [ Z W((L, ) )(P.)/P] (12) 

A=1 A ^ 


IV. 2 Constraints 


Airport authorities seeking to reduce noise and still service passenger 
demand for their airports may impose various related operating constraints. 
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IV . 2 ♦ 1 FHfiht Limitations 

For any given airport flight limitations may be imposed for many dif^- 
ferent reasons and consist of many different constraints on aircraft opera- 
tions. The flight limitations of primary Interest in this research are 
those operational constraints that may be imposed in on attempt to improve 
the noise environment around an existing airport. The operational flight 
constraints formulated in this research are: 

1. The annoyance as measured by a specified community response 
index may not exceed critical values for specified communities 
in the airport vicinity, 

2. The number of certain types of aircraft that may operate into 
and out of a given airport cannot exceed some upper limit, 

3. The number of certain types of operations such as the take-offs 
of long stage length aircraft may be limited. 

4. Certain runways may only be used during certain periods or a 
limitation on the number of operations per period for any 
given runway may be specified. 

5. Certain trajectories that correspond to specified ground tracks 
may only be used during certain periods or a limitation on the 
number of operations per period for any given ground track may 
be specified, 

6. The total number of operations per time period may be constrained. 

These flight limitations are now formulated mathematically: 

1. The annoyance as measured by some community response index 
for a given area, A, must not exceed some critical value, 


RI, 


i ''a 


A - 1 


> » ? • J 


NA 


(13) 
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For example if was selected as the community response index 
then (13) would appear as a function of the decision variables 
as 


NI NJ NK (NEL,,,./10) 

10 log,. S S S {(X,,,, +10X,,,„)10 ^ 

i«l j«l k«l 


+ (v^ + lOYy^JlO 


(NEL,,,/10) 

} - 49.4 ■: V, 
- A 


A « 1, , . , , NA 


(14) 


2. The number of aircraft type t allowed to operate in and out of a 
given airport is limited to N operations per day 

U 


NK NJ NL 

i: z z X. 


NJ NL 


‘tki£ ^ ^ ^ti; - ^t 

k=l j=*l A=1 j=l SL^l 


t = 1, ..., Nl (15) 


3. The number of take-offs with stage length s is limited to N^ 
operations per day 


NI 

NJ NL 




Z 

1=1 

Z Z 
j=l il=l 

X., < N 

IsjJl - s 

s = 1, NK 

(16) 

The 

number of 

operations for 

runway r is limited 

to N^ per day 

NI 

NK 

NL 

NI NL 


Z 

i=l 

Z Z 

k=l jsR. 
J 


Z Z Z Y < 

i=l jeR, i=l " 

N , 
r 


r = 1 


> • • • > 


NR 


(17) 
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5. The number of operations corresponding to ground track d is 
limited to per day 


NI NL 

t £ Z X + Z Z Y, < N, d ° 1. NJ (18) 

1 k A 1«1 JI--1 ^ 


6. The total number of operations in time period p must not exceed 
N 


NI 

Z 


NK 

Z 


NJ 

Z 


i“l k^l j'=l 


X + 
ikjp ^ 


NI 

Z 


NJ 

Z 


i°l j<=l 


ijP 


< N 


1 , 


NI 


(19) 


IV, 2, 2 Aircraft Availability 

Only a limited number of the various types of aircraft servicing an 
airport will be available during each time period for either landing or 
departure. This may be expressed analytically as 

NJ 

Z < NY^^ i = 1, NI; f <= 1, NL (20) 

j=l 

m NJ 

Z 2 X,. < NX 1 ^ 1, NI; f = 1, NL (21) 

k=l j=l 

IV, 2. 3 Passenger - Aircraft Demand 

Passenger demand may be established for any given airport. Then 
passenger demand may be translated into aircraft demand. Such demands may 
take the form of requiring at least N^^ operations of interest along a sub- 
set of tracks (where n denotes the subset of interest) during time period H. 
For departures this may be represented by, 
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NI NK 

K S X X,,.. > N, £ “ 1, NL; n « 1, NN (22) 

i«I k=l ^ 

where, NN ° the number of track subsets. 

This is one of the simplest ways in which passenger demand may be accounted 
for. More eleborate demand relationships could be derived and utilized if 
desired, 

Collectively Equations (12) through (22) define a mathematical model 
for this research. This model may be classified as a nonlinear integer 
mathematical programming model. Solution techniques for such a model are 
discussed in the succeeding section. 

V. SOLUTION TECHNIQUES AND COMPUTATIONAL RESULTS 

Examination of the mathematical model formulation in Section IV reveals 
that the objective function is nonlinear. The constraints, with the excep- 
tion of (13), are linear. It may be classified as a nonconvex programming 
problem of from 30 to 200 constraints with 100 to 500 variables for small to 
medium-small airports. Only the smallest of such problems would even be 
attempted through the direct application of one of the existing nonlinear 
optimization algorithms. Even if the attempt were made there would be no 
guarantee that the global optimum would be identified. The applications of 
interest may give rise to mathematical models with as many as a few thou- 
sand variables and several hundred constraints. The only optimization solu- 
tion techniques that appear to be feasible for application to such size 
problems would require linearization (approximation) of all nonlinear 
equations. 
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V . 1 Selected Solution Alsorlthm 

The ooluUion technique selected by this author is to use successive 
linear approximation for the nonlineor objective function and to consider 
the decision variables ns continous. This algorithm, summarized below, 
uses gradient information to form successive linear objective function 
approximations. The objective equation in the following Step 1 is a 
linear approximation to the sum of all the population weighted L^^values 
ar- ii function of the control variables. 

Step 1 

^ NA 

minimize Z « I (23) 

A°1 ^ ^ 


subject to all linear constraints 


S 


A 


a 


NI NK NJ 
SEE 
i«l k«l j=l 


NEL 


ikjA _ 


{10 


10 


4.94) 

^^ikjl ^°^ikj2^ 


NEL, , . 

- 4.94) 

+ 10 


(24) 


letting, 


‘ikjA 


10 


(JA. 4.5,) 


NEL 


ijA _ 


and 


ijA 


10 


10 


4.94) 


NI NK NJ 
SEE 
i=l k=l j=l 


^^ikjA^^ikjl 


lOX,. .„) 
ikj 2 ' 


implies 
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From Equations (1^) and (25) we see that 




( 25 ) 


(26) 


which means that the set of nonlinear constraints in (14) may be handled by 
introducing the additional set of variables S ^^^and upper 

bounding them, i.e., 

(V./IO) 

< 10 A “ 1, . . ., NA (27) 


This will provide a feasible solution to the problem but in no sense 
guarantees an optimum. Instead of the function given in Equation (23), the 
nonlinear objective function should be derived from Equation (12) upon sub- 
stituting Equation (10) and Equation (26) . 


minimize 


,N 


NA 3.364 X 

^ O ^"-^TTS 

A=1 + 1.43 X 10 


(28) 


However the Step 1 objective function has provided very good solutions to 
the original problem for the few example problems solved. 

The reason for such an objective function is that it is linear in the 
decision variables and subject to linear constraints, hence corresponds to 
a linear programmiitg problem. With sophisticated computer implementation, 
linear progranm:ii<g solution techniques are capable of solving very large 


problems, For example, problems Involving several thousand linear constraint 
equations and tens of thousands of variables are within the realm of possibility 
for the very efficient linear programming solution algorithms implemented on 
modern computers. However it should also be pointed out that to produce 
solutions to such large problems requires extensive efforts in data preparation, 
manipulation of the linear programming computer code, and interpretation of 
computational results, 

Step 2 

Obtain a truncated Taylor series expansion about the solution point 
from Step 1, say S_* 

u(^) = Z^(S*) + (S - Sf<) (29) 

Now minimize (29) subject to the original linear constraints. This 

corresponds to solving another linear programming problem. Denote the 

1 

solution to this lliiear program s29) as ^ . Since u(^) is constructed from 
N 

the gradient of Z at an improved solution point can be secured only 

if u(^^) < u(£*). This will not guarantee that Z^(£^) < Z^(£*) unless ^ 

1 

is in the immediate neighborhood of _S*. However, given u(^ ) < u(^*) there 

2 1 
must exist a point, say ^ , on the line segment between ^ and _S such 

that Z^(^^) < Z^(£*). To determine ^ one solves 

minimize Z^^(^* + a(^^ ^*)) 

a 

Z^(S^) = Z^(S* + - S*)) = minimize Z^(_S* + a(^^ S*)) (30) 

0<a<l 


2 

Set S* equal to _S and repeat Step 2 as many times as required to obtain 

1 

the stopping condition, u(S ) > u(&*). At this point no further improve^- 
ment is possible. 

The linear programs generated at the successive iterations of this 
algorithm differ only in the coefficients of the objective function. There- 
fore, the sensitivity analysis options available on a computerized linear 
programming solution algorithm mav be exercised to efficiently carry out 
all calculations required for the second through the last iteration. The 
solution to (30) may be obtained by the use of any one-dimensional search 
technique. Reference [9] provides a description of several one-dimensional 
search algorithms that could be used. 

The MPS III (Mathematical Programming System) was utilized to imple- 
ment this solution procedure on an IBM 370 computer. DATAFORM, a data 
managevrtent subsystem available through the MPS III package, was used to 
generate data, Interface the required FORTRAN programs and produce computa- 
tional reports for all application airports considered. The next section 
of this paper describes the application of the implemented solution pro- 
cedure to an example airport, 

V. 2 Application Airport for Solution Procedure 

The following describes a medium sized airport to which the solution 
procedure, just described, has been applied. The analysis was limited- to 
commercial airline traffic for the particular airport of interest, A series 
of tables and figures, presented in Appendix A, are used in constructing the 
appropriate mathematical model. The computational results obtained from 
application of the solution algorithm to this example airport are displayed 
in Tables 5 and 6. Noise minimal operating scenarios for several other air- 
ports have been computed using this same solution algorithm. However, 
because of length considerations they are not presented in this paper. 
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TABLE 5 

COMPUTATIONAL RESULTS FOR APPLICATION AI1U?0RT (DEPARTURES) 
0.17723 


Number of A lrcraf L' Stage Track T1 

Operations DC-9 [*727 [ DC-S T L-1011 Length Number Da 


COMPUTATIONAL RESULTS FOR APPLICATION AIRPORT (ARRIVALS) 





^ St* s<^ KA Sr* s>* S^ Sd* s«* 












0.14457 


Number of 
Operations 


TABLE 6 

COMPUTATIONAL RESULTS - Example 2 

Departures 


Track 

Lensth I Number 


Aircraft 

DC- 9 

727 

X 


X 

X 


X 

X 


X 


X 

X 


X 


TABLE 6 (Cont'd) 

COMPUTATIONAL RESULTS (ARRIVALS) 


X 


X 

X 


X 

X 


X 


X 

X 


X 


5 

X 

5 


6 i 

• X 

6 

i 

7 

X 

7 


8 

X 

8 

X 

8 

i 

i 


1 ^, 
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TABLE 8 

COMPUTATIONAL RESULTS - Example 4 
Nil 0.27505 Departures 


Number of 
Operations 

Aircraft 

HSU 


Time Period 

DC- 9 

727 

A- 300 

DC- 10 


mmm 

71 

X 




1 

1 

X 


44 

X 




2 

1 

X 


23 


X 



3 

1 

X 


2 


X 



1 

1 


X 

1 


X 



2 

1 


X 

1 


X 



3 

1 


X 

71 


X 



1 

2 

X 


44 


X 



2 

2 

X 


23 


X 



3 

2 

X 


2 


X 



1 

2 


X 

1 


X 



2 

2 


X 

1 


X 



3 

2 


X 

7 

X 




1 

3 

X 


64 


X 



1 

3 

X 


44 


X 



2 

3 

X 


23 


X 



3 

3 

X 


2 


X 



1 

3 


X 

1 


X 



2 

3 


X 

1 


X 



3 

3 


X 

71 


X 



1 


X 


24 


X 



2 


X 


11 


X 



3 


X 


20 



X 


2 


X 


12 




X 

3 


X 


2 


X 



1 



X 

1 


X 



3 



X 

1 



X 


2 



X 
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TABLE 8 (Cont’d) 

COMPUTATIONAL RESULTS - Arrivals 



















TABLE 9 

COMPUTATIONAL RESULTS - Example 5 
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TABLE 10 

Annoyance Comparison for Application Airports 



Minimum 

Operating 

Conditions 

Current 

Operating 

Conditions 

Example 1 



Noise Impact Index (Nil) 

0.17723 

0,25648 

^Number of Persona Highly 
Annoyed in 24 Ibur Period 

36,578 

52,937 

Example 2 



Noise Impact Index (Nil) 

0.14457 

0.16792 

^Number of Persons Highly 
Annoyed 

4775 

5546 

Example 3 



Noise Impact Index (Nil) 

0.15479 

0.18520 

^Number of Persons Highly 
Annoyed 

39,352 

47,082 

Example 4 



Noise Impact Index (Nil) 

0.27505 

0.31086 

Number of Persons Highly 
Annoyed 

744,385 

841,286 

Example 5 



Noise Impact Index (Nil) 

0.28658 

0.37177 

^Number of Persons Highly 
Annoyed 

126,250 

163,779 


^These quantities were calculated using the formula on 
page B“5 of reference [15]. 
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TABLE 6 


Annoyance Comparison for Application Problem 



Minimum 

Current 


Operating 

Operating 


Conditions 

Conditions 

Noise Impact Index (Nil) 

0.17723 

0.25648 

* Number of Persons Highly 
Annoyed in 24 Hour Period 

36,578 

52,937 


’^These quantities were calculated using the formula on page 
B-S of reference [15]. 


One can note from Table 6 that an estimated 31 percent reduction in 
the number of people highly annoyed may be achieved by utilizing the noise 
mlnlmlal operating procedure, This reduction in Impact is typical of all 
airports analyzed, l.e., reductions of approximately 20-40% are achievable 
at the airports analyzed to date. 

VI . CONCLUSIONS 

An optimization mathematical model whose objective is to minimize a 
measure of annoyance due to the arriving and departing aircraft for a given 
airport has been formulated. A corresponding solution algorithm, relying 
upon the solution of linear programming problems, was subsequently developed 
and computational results for one operational commercial airport was pre- 
sented, The solution algorithm, even though it does not guarantee to find 
the global optimum, should produce very good solutions for any given airport. 
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Our experience has been a 10-40 percent reductions in noise impacts ore 
possible. This includes several airports that hod previously implemented 
a noise preferential runway assignment system. Such systems are designed 
by airport authorities to reduce noise exposure taking into consideration 
population distributions, existing geographical features that might 
provide natural corridors (e.g., a river), airport Instrumentation, 
weather conditions, etc. An especially attractive feature of the 
suggested solution algorithm is that it is capable of solving very large 
problems. For example, it would be feasible to attempt the solution of 
problems involving several thousand variables and 500 plus linear con- 


straints. 
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APPENDIX A 

Data for Appllcatloi* Airport 


TABLE A.l 

Types of Aircraft Considered 


Aircraft Type 

Number of Stage Lengths* 

DC-9-32 

3 

B727-200 

4 

DC-8-55 

4 

L-1011 

4 


* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles) , 2 (500-1000) nautical miles) , 3 (1000 
-1500 nautical miles) and 4 (1500-2500) 


TABLE A. 2 

Demand for Incoming Flights 


Runway Tracks 

Day 

Night 

(15)24* 

6 

1 

(13)12 R 

90 

10 

(14)30 L 

98 

11 


* The number in parenthesis designates the track number de- 
fined in Figure A.l. The following number designates the 
associated runway. 



Avail a 


■pe. 
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FIGURE A. 2 3° Approach Profile. ("Airport Noise Reduction 

Forecast Volume II — NEF Computer Program Descrip- 
tion and User's Manual." Study for Office of Noise 
Abatement, Department of Transportation by Wyle 
Research. Washington, D.C., October 1974.) 


C 



FIGURE A. 3 Take-off Profile, ("Airport Noise Reduction Fore- 
cast Volume II — NEF Computer Program Description and 
User's Manual." Study for Office of Noise Abatement, 
Department of Transportation by Wyle Research. 
Washington, D.C., October 1974.) 
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FIGURE A. A 

Population Areas in the Vicinity of Example Airport 
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TABLE A, 6 

Population of Areas in Vicinity of Example Airport 


Area 

Pop 

Area 

Pop 

1 

20078 

34 

879 

2 

16176 

35 

2127 

3 

8022 

36 

1161 

4 

8732 

37 

2591 

5 

11887 

38 

7447 

6 

12317 

39 

4294 

7 

4987 

40 

3585 

8 

5822 

41 

4401 

9 

7053 

42 

8060 

10 

19680 

43 

6075 

11 

9494 

44 

13940 

12 

3579 

45 

9578 

13 

5596 

46 

10334 

14 

946 

47 

6659 

15 

96 

48 

17991 

16 

8918 

49 

17268 

17 

5339 

50 

6714 

18 

9475 

51 

3048 

19 

532 

52 

8144 

20 

8874 

53 

13093 

21 

3799 

54 

5193 

22 

3996 

55 

5359 

23 

453 

56 

21192 

24 

1250 

57 

13785 

25 

1570 

58 

5640 

26 

4109 

59 

16827 

27 

5833 

60 

17408 

28 

6908 

61 

7977 

29 

107J6 

62 

15239 

30 

11813 

63 

10034 

31 

2399 

64 

36311 

32 

15610 

65 

10852 

33 

10661 
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APPENDIX B 

Data for Application Airport 
Example 2 


TABLE B.l 

Types of Aircraft Considered 


Aircraft Type 

_ ... 

Number of Stage lengths* 

DC-9-32 

1 

B727-200 

1 


* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
“1500 nautical miles) and 4 (1500-2500) 


TABLE B.2 

Demand for Incoming Flights 


Runway Tracks 

Day 

Night 

(5) 03* 

5 

1 

(6) 21 

5 

1 

(7) 15 

5 

1 

(8) 33 

5 

1 


* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
-1500 nautical miles) and 4 (1500-2500) 
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APPENDIX C 

Data Application Airport 
Example 3 

TABLE G.l 

Types of Aircraft Considered 


Aircraft Type 

Number of Stages* 

DC-9-32 

2 

B727-200 

2 

A-300 

1 


* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
-1500 nautical miles) and 4 (1500-2500) 


TABLE C.2 

Demand for Incoming Flights 


Runway Tracks 

Day 

Night 

(5) 12R* 

33 

4 

(6) 30L 

33 

4 

(7) 03R 

33 

4 

(8) 21L 

33 

4 


* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
-1500 nautical miles) and 4 (1500-2500) 



TABLE C,3 

Take-Off Flight Demands 


Tracks 

Stage 

Stage 1 


1 

2 

1 

2 

(1) 12R 

23 

9 

2 

1 

^2) 3QL 

23 

9 

2 

1 

(3) 03R 

23 

9 

2 

1 


23 

9 

2 

1 


TABLE C.4 

Available Aircraft for Arrivals 


DC-9-32 

B727-100 

-300 


48 

2 

84 

15 

1 

.1 


DG-9-32 

B727-100 


TABLE C.5 

Available Aircraft for Departures 
a 


48 

2 

84 

15 

1 

1 




















































































































APPENDIX D 

Data for Application Airport 
Example 4 


TABLE D.l 

Types of Aircraft Considered 


Aircraft Type 


Number of Stages* 


DC“9-32 

B727-100 

A-300 

DC-10-10 


* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
-1500 nautical miles) and 4 (1500-2500) 



TABLE D.2 

Demand for Incoming Plights 


Night 


139 

6 

139 

6 

139 

6 

139 

6 




* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
-1500 nautical miles) and 4 (1500-2500) 
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APPENDIX E 

Data for Application Airport 
Example 5 
TABLE E.l 

Types of Aircraft Considered 


Aircraft Type 

Number of Stages* 

DC-9-32 

3 

B727-100 

4 

B707-320 

6 

A-300 

3 

L-1011 

6 

DC-10-10 

6 


* Tlie stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
-1500 nautical miles) and 4 (1500-2500) 

TABLE E.2 

Demand for Incoming Flights 


Runway Tracks 

Day 

Night 

(7) 9R* 

96 

17 

(8) 27L 

96 

17 

(9) 9L 

96 

17 

(10) 27R 

96 

17 

(11) 12 

96 

17 

(12) 30 

96 

17 


* The stage length of a departing aircraft is a measure of 
the distance to the next destination. Stage lengths 1 (0- 
500 nautical miles), 2 (500-1000) nautical miles), 3 (1000 
-1500 nautical miles) and 4 (1500-250?.*; 
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TABLE E.3 




Take-Off Flight: Demands 


Runway 

Tracks 


Day 

Night 


Stage 

Stage 


1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

(1) 

9R 

29 

23 

35 

5 

1 

1 

6 

1 

4 

2 

1 

1 

(2) 

27L 

29 

23 

35 

5 

1 

1 

6 

1 

4 

2 

1 

1 

C3) 

9L 

29 

23 

35 

5 

1 

1 

6 

1 

4 

2 

1 

1 

(4) 

27R 

29 

23 

35 

5 

1 

1 

6 

1 

4 

2 

1 

1 

(5) 

12 

29 

23 

35 

5 

1 

1 

6 

1 

4 

2 

0 

1 

(6) 

30 

29 

23 

35 

5 

1 

1 

6 

1 

4 

2 

0 

1 


TABLE E.4 

Available Aircraft for Arrivals 


Type 

Day 

Night 

DC-9-32 

155 

15 

B727-100 

271 

42 

707-320 

47 

23 

A-300 

12 

1 

L-1011 

16 

13 

747-100 

17 

9 


TABLE E.5 


Aval] 

Lable Aircraft for Departures 

Type 

Day 

Night 

DC-9-32 

155 

15 

B727-100 

276 

42 

707-320 

47 

23 

A-300 

12 

1 

L-1011 

16 

13 

747-100 

17 

9 
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